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Abstract 

> ; 

It is very plausible that confinement in QCD comes from linking of Wilson 
O loops to finite-thickness vortices with magnetic fluxes corresponding to the 

center of the gauge group. The vortices are solitons of a gauge-invariant QCD 
00 ■ action representing the generation of gluon mass. There are a number of other 

solitonic states of this action. We discuss here what we call nexus solitons, in 
J . which for gauge group SU(N), up to N vortices meet at a center, or nexus, 

provided that the total flux of the vortices adds to zero (mod N). There are 
qj \ fundamentally two kinds of nexuses: Quasi- Abelian, which can be described 

^ ■ as composites of Abelian imbedded monopoles, whose Dirac strings are can- 

celled by the flux condition; and fully non-Abelian, resembling a deformed 
^ . sphaleron. Analytic solutions are available for the quasi-Abelian case, and 

we discuss variational estimates of the action of the fully non-Abelian nexus 
solitons in 577(2). The non-Abelian nexuses carry Chern-Simons number (or 
topological charge in four dimensions). Their presence does not change the 
fundamentals of confinement in the center- vortex picture, but they may lead 
to a modified picture of the QCD vacuum. 
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I. INTRODUCTION 



The center-vortex picture of confinement in QCD was proposed long ago 
jl], 0, |3], f|, |]. In essence it identifies the QCD vacuum as a condensate of 
finite-thickness vortices of co-dimension 2 in d=3,4; these vortices (closed 
loops in d=3, closed surfaces in d=4) have a finite transverse extension mea- 
sured by the inverse gluon mass and possess magnetic fluxes which are ex- 
pressed via elements of the center of the gauge group. 

In more recent times there has been a resurgence of work on the center- 
vortex picture, primarily on the lattice || || |Hj. This has resulted in 
spectacular confirmation of the center- vortex picture of confinement, with 
the SU(2) string tension being fully captured by by the vortex-dominance 
approximation of replacing the exact Wilson loop value by (— ) N , where N 
is the total linking number of vortices linked to the loop. 

The center vortices resemble Dirac strings in carrying quantized magnetic 
flux confined to a closed loop or surface, but differ radically from them in 
that the center-vortex flux is spread out over a finite transverse extent. If the 
flux were instead fully localized (to a delta-function in the continuum, or to 
a single lattice spacing on the lattice) two things would go wrong: First, the 
action of the vortices would be quadratically infinite, either in the continuum 
[p]] or on the lattice . Second, these thin vortices would contribute nothing 
whatsoever to the adjoint "string tension" (more properly, a perimeter term 
in the potential resembling a linearly- rising potential which breaks), leaving 
only a Coulomb term. In fact, lattice simulations [jllj , |i~2l and continuum 
theory O, [14], [L5[ show that the thick vortices of the center- vortex picture do 



lead to a breakable "string" potential for the adjoint representation, directly 
attributable to the thickness of the vortices which can only partially overlap 
the Wilson loop. 

Avoidance of Dirac-string singularities is a key point in constructing soli- 
tons in the continuum center-vortex picture, and will be important for us 
in constructing the objects considered here, which we call nexus solitons. A 
nexus soliton has several thick vortices (up to N of them for gauge group 
SU(N)) which meet at a central region which is essentially of magnetic- 
monopole character, with exponentially-falling field strengths as one moves 
out from the center. In order that this central region have finite action, the 
magnetic fluxes of the several vortices meeting at the center should add to 
zero. If not, the monopole region has fields which diverge like r~ 2 at the ori- 
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gin, like the Wu-Yang monopole or the broken-symmetry monopoles found 



by Corrigan et al. long ago ||16|| . Moreover, if the fluxes do not add to zero 
there can be long-range monopole fields behaving like r~ 2 at large distances, 
as well as naked Dirac strings. These last two phenomena are due to the 
termination of Dirac strings at the monopole, in the case where the fluxes do 
not add to zero (mod N). 

In previous solitons of the center-vortex picture Jl| this unwanted be- 
havior is simply avoided, in part by having the loop or surface representing 
the vortex be closed, and in part by cancellation of Dirac-string singulari- 
ties between a transversely-extended part of the solution and a naked Dirac 
string. (In fact, the center vortex is simply an Abelian Nielsen-Olesen vortex 
embedded in the gauge group but with infinite Higgs mass, or otherwise de- 
scribed it is a soliton of an effective low-energy QCD action consisting of the 
usual G 2 term plus a gauged non-linear sigma model term, giving a gauge- 
invariant mass to the gluon.) Any boundaries of the loop or surface gives 
rise, as is well-known, to objects mathematically identical to the original 
Dirac monopole plus associated string. In addition to these closed vortices, 
the center-vortex picture also has a spherically-symmetric static solution, re- 



sembling a sphaleron but with no symmetry breaking [17]. It has no Dirac 
strings in the usual spherical gauge presentation.^ 



The existence of nexus solitons was alluded to long ago [1], [L8| , but to the 
author's knowledge these have not been explicitly constructed. We partially 
fill that gap here, both for SU(2) and SU(3). In certain cases, which we call 
quasi- Abelian, we can give explicit analytic results and in other cases, which 
are fully non-Abelian, we restrict ourselves to variational estimates. These 
nexus solitons go considerably beyond the original vortex solutions, which 
as mentioned above are really imbedded Abelian Nielsen-Olesen vortices; 
the new solitons really cannot be described by Abelian embeddings. The 
non-Abelian solution which we discuss is essentially a QCD sphaleron |fPT| 
from which several fat vortices emerge. It carries topological charge (or 
Chern-Simons (CS) number in d=3). It can be thought of as a sphaleron- 
vortex composite, and in SU{2) is energetically stable against decay into 
a sphaleron and a vortex. In SU(3) the quasi-Abelian nexus with three 
unit-flux vortices is stable against decay into a QCD sphaleron on energetic 

^ike the Wu-Yang monopole it can be transformed into an Abelian gauge where it 
does have a Dirac string, but the required gauge transformation is singular. 
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grounds; this stability is not likely for the SU(3) non-Abelian nexus, which 
is much heavier. 



It is worth noting [H5 that even the original Abelian center vortices can 
carry topological charge, which in d=3 is associated with the linking or self- 
linking of vortices (in d=4, this is expressed as an intersection number of 
surfaces). This topological charge is not necessarily integral, but can be 
quantized in units of 1/N for gauge group SU(N). This need not be a con- 
tradiction with periodicity of 2% in the 9 angle, or with integral quantization 
of the Chern-Simons level. 

Aside from questions of ^-dependence of QCD dynamics, which we will 
not discuss here, of what good are these vortices? The answer is, for confine- 
ment at least, that they are not good for much of anything; they do not [|TJ 
change the basic picture of confinement for the fundamental representation 
which arises when one uses center vortices without nexus solitons. An inter- 
esting question which we are not in a position to answer yet is whether nexus 
solitons, which are capable in principle of internetting the center-vortex vac- 
uum in a very complicated way (see Fig. 4(b) below), actually do so. Put 
another way, the question is whether the change in entropy arising from this 
internetting is in a sense and of a size which can over come the positive action 
needed for the nexus monopole regions. If this happens the vacuum conden- 
sate will look quite different from the usual (naive) center- vortex condensate, 
and possibly could be detected in lattice calculations. We will discuss in the 
concluding section what could be learned about nexus solitons and other 
elements of the center-vortex picture from the lattice. 

We begin with a quick review of the elements of the center- vortex picture, 
including the effective action. 



II. THE CENTER-VORTEX PICTURE 



For convenience we will from now on, except when otherwise noted, work in 
three dimensions, since all the solitons we find are static. 

It was argued long ago [T7| that the low-energy effective energy functional 
E of QCD consisted of the usual three-dimensional action plus a gauged 
non-linear sigma model term, representing the generation of a gluon mass. 
In turn, the gluon mass appears pO] because of the fundamental infrared 



instability of QCD in d=3,4, where the Schwinger- Dyson equations only have 
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non-singular solutions when a mass of sufficient size is generated. The energy 
isf] 

E = -j J d 3 x{^-Tr[G i:i (x)] 2 + M 2 Tr\D t U\ 2 }. (f ) 

This action is usable only at low energies, because at high energies we must 
account for the vanishing [^0], |21] of the mass, something like M 2 (q) ~ 
(G 2 )/q 2 where (G 2 ) is the condensate expectation value.f] The mass must 
vanish at large momentum if we are to find solutions to the Schwinger-Dyson 



equations 20 



The fundamental Abelian- vortex solution of the equations of motion com- 
ing from the equations of motion of (1) is: 

Ai(x; J) = (2nQ J /i)e ijk d j I dz k {A M (x - z) - A (x - z)} (2) 

Here the integral is over a closed loop of coordinates z,Am (A )is the massive 
(massless) free propagator in three dimensions, and Qj is an N x N SU(N) 
matrix in the Cartan subalgebra. It is normalized so that the long-range 
(A ) part of the Wilson loop integral P exp / dxi-Ak(x) for a Wilson loop 
linked once to the vortex is in the center of the group. For the fundamental- 
representation Wilson loop, this means 

exp(2vriQj) = exp(2vri J/N). (3) 

A convenient set of matrices Qj is formed from a set of fundamental Cartan 
matrices Qf. 

Q t = diag(l/N, . . . 1/N, -1 + 1/N, 1/N, . . . 1/N) (4) 

where — 1 + 1/ N is in the zth position. Only N—l of these matrices are linearly 
independent, because the sum of all of them is zero. Then the matrices 



2 We use the usual anti-hermitean gauge potential matrix (A a /2i)A a (x), with Tr\ a X b = 
(l/2)6 a b', Di = di + Ai. The matrix U is a unitary matrix undergoing left transformations 
when the gauge potential is gauge-transformed. 

3 We will not actually account for this high-energy decrease in mass, with the conse- 
quence that the second term in the action (1) is logarithmically divergent at short dis- 
tances, for the vortices. Of course, if the mass vanishes quadratically at short distances in 
actuality this is a spurious divergence. In fact, the nexuses themselves do not suffer from 
this logarithmic divergence. 
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Qi + Qj + Qk + ■ ■ •, with i 7^ j 7^ k ^ ■ ■ -, obey (3) if there are J terms in the 
sum. For 1 < J < [N/2], where [•] indicates the integral part, we can choose 
in any convenient way one representative of Qj from such sums, and for 
greater values of N we use instead the matrices Qn-j = Q-j, representing 
anti- vortices of flux —7. For future use we record: 



)2 _J(N-J) 



Tr Qj = v N ' (5) 

Confinement now arises in d=3,4 [[]J by averaging over the fluctuating 
phase factors, equation (3), which arise in evaluating the expectation value 
of the Wilson loop. 

For ££7(2) there is only one non-trivial vortex, and it is self-conjugate; 
for £77(3) there is one vortex and its anti-vortex. 



III. NEXUS SOLITONS 

We are looking for solutions to the equations of motion coming from (1) 
which have several thick vortices meeting at a monopole-like center, or nexus, 
with finite action and non-singular field strengths. For ££7(2) the only non- 
trivial case is when two J = 1 vortices meet, as in Fig. 1. By convention we 
choose ingoing arrows, representing vortex magnetic fields directed toward^ 
the nexus, to describe a nexus, and the opposite direction corresponds to 
an anti-nexus. For £77(3) there is only one additional non-trivial nexus, as 
given in Fig. 2, in which three 7=1 vortices meet. 

Actually, these figures are slightly misleading, since ultimately the vortex 
strings must close. The true situation, therefore, is that of Figs. 3 (for ££7(2)) 
and 4 (for ££7(3)) and generalizations thereof, which have as many nexuses 
as anti-nexuses. We note a device for constructing pictures of nexuses: the 
££7(3) graphs of Figs. 2 and 4 are topologically the same as baryonic Wilson 
loops, with baryon violation of ±1 unit at every nexus or anti-nexus (given 
that the baryon number of a quark is 1/3). Of course, such B- violation 
is entirely absent in QCD, which is vector-like, and this is only a heuristic 
analog. 

4 Defining the sense of magnetic fields in non-Abelian theories requires some conven- 
tions. We choose to assign a positive sign to all the Qj and a negative sign for Qn-j = Qj 
for 1 < J < [N/2]. 
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We begin by looking for nexuses which are essentially Abelian, like the 
center vortices themselves. There are two possible types: U(l) nexuses, 
and quasi- Abelian SU(N) nexuses with charges in the Cartan subalgebra. 
In both cases the only possible nexuses are constructed of thick vortices 
whose fluxes add to zero at the nexus core. However, the U(l) nexuses are 
marginally stable to decay into simple loops, which is not so for the quasi- 
Abelian nexuses. 

It will be evident that the essentially Abelian nexuses we find in this sec- 
tion are in some sense composites of monopole-like excitations. However, it 
will turn out to be impossible to separate the nexuses into their fundamental 
parts, since each part is a monopole with a thick vortex carrying the flux 
into the central monopole region, but a naked Dirac string is needed to carry 
the flux out. In a nexus all these naked Dirac strings cancel because of the 
condition that the total flux is zero; they do not cancel for the pieces from 
which the nexus is made. 



A. U(l) nexuses? 

We ask whether there is a nexus in a U(l) gauge group, like electromagnetism, 
but with a mass term; the physical analog would be a superconductor. It 
should be apparent that all the vortex fluxes must add to zero at the nexus, or 
else there will be monopole fields which are either singular at short distances 
or are long-range or both. However, for future use it is worth detailing the 
explicit forms of the difficulties encountered with non-zero total flux. We first 
give the well-known form of the solution (2) for an Abelian center vortex lying 
along the z-axis, in standard cylindrical-coordinate notation: 

A l (x) = ^ t [MK l (Mp)--]. (6) 
i p 

Here K\(Mp) is the Hankel function of imaginary argument, which falls off 
exponentially at large argument but behaves like 1/Mp near p — 0. This term 
comes from the Am propagator in (2), while the 1/p term comes from A ; 
this latter term is a long-range pure-gauge term, obviously singular at p = 0. 
The K\ term just cancels the explicit 1/p term. This is what we meant earlier 
when we said that the center vortex has an extended part with a Dirac string 
just cancelling an explicit Dirac string. If this did not happen there would 
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be a field strength proportional to S(x)S(y), which yields a quadratically- 
divergent action. Clearly this cancellation occurs for any closed-loop integral 
in (2), since A M behaves like A at short distances. 

Now consider terminating the string integral in (2) at some convenient 
point, say the origin. Then the A^ part of (2) gives a thick vortex along the 
negative z-axis plus a massive monopole whose (isotropic) field strength goes 
like 1/r 2 at the origin but decays exponentially at infinity. The A part gives 
the original Dirac monopole, with its string cancelling that of the massive 
term as before, but also with a long-range monopole with field strength ~ 
1/r 2 at both short and long distances. The short- distance part of this field 
strength is cancelled by that of the massive term, but the long-distance part 
survives. This does not happen in QCD, where all field strengths (but not 
all gauge potentials!) are short-ranged, so it would seem that the A part of 
(2) must have an integral over a closed loop. But then if we terminate the 
massive part of (2) at some point, its short-distance field behavior of 1/r 2 
is uncompensated, and leads to a badly-divergent action. Moreover, there is 
now a naked Dirac string coming from the long-range part of ((2). All of this 
is unacceptable. 

Of course, the way to cure these problems is to make sure that the sum of 
vortex fluxes is zero, whenever several strings terminate at at common point. 
The gauge potential is of the form 



A{x) = —tijkdj YJlA / dz k (A M (x - z) - A (x - z))\ (7) 



where the sum of the U (1) fluxes qA is zero, and the Ta are string paths which 
meet at the origin, as in Fig. 5. Each path integral runs from to infinity. 
All the long-range monopole fields from the A terms cancel, because they 
are isotropic and string-independent; the short-range parts of these terms 
cancels the 1/r 2 field strengths from the massive terms at the origin, for the 
same reason; and the Dirac strings all cancel between the short-range and 
long-range terms, just as for the center vortex. 

Note that this is not useful as a model for SU (2) nexuses, since with two 
fluxes q\ = —q 2 the fields do not correspond to Fig. 1, but rather one field 
points into and the other out of the nexus. This just gives the usual center 
vortex with no nexus, at least if the strings joining at the nexus lie along 
a common axis, say the z-axis. The SU(2) nexus is truly non-Abelian, and 
will be discussed in Section IV. 



i 
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There is another more useful interpretation of the gauge potential (7), 
in which the strings form closed loops, some segments of which lie along a 
common line. Examples of triple nexuses, corresponding to SU(3), are given 
in Fig. 6, which may be considered as possible realizations of Fig. 4. In Fig. 
6(a), if the strings labeled 1,2 have incoming fluxes qi,q2 respectively, which 
we will take to be positive, then the double string has incoming flux —qi — q 2 
and automatically respects the condition of zero flux sum. However, the 
configuration of Fig. 6(a) is manifestly unstable into decay into its consituent 
loops, because the action of the double string, proportional to (q± + q 2 ) 2 , is 
larger than the sum of the actions for strings 1 and 2. 

In Fig. 6(b), if the string fluxes are chosen to add to zero, then the triple 
string down the middle is a phantom string, with no energy, because its flux 
is zero. This configuration is not unstable to decay into its constituent loops, 
which would require supplying the energy needed to separate the phantom 
loop into three real loops. However, it is not a useful description of true 
SU(3) nexuses, because these will have the same action (per unit length) 
along each of its three vortices. In a U(l) group where the qi are three real 
numbers, it is impossible to find three real non-zero numbers whose sum is 
zero and whose squares are all equal. 

We conclude that a U(l) description of nexuses is not possible, and turn 
to quasi-Abelian nexuses. 

B. Quasi-Abelian nexuses 

The only interesting SU (2) case is the non-Abelian case of the next sec- 
tion, so we start with £77(3). We use the form (7) 

27r 3 r 

Mx) = —Y / [T aa dz k (A M (x - z) - A (x - z))} (8) 
1 0=1 ■ /r ( a ) 

where T n , T 22 , T 33 are defined in the Appendix; these matrices, whose sum is 
zero, are just the flux matrices Qi of equation (4) for £77 (3) but with different 
names. They obey 

exp(27iiT aa ) = exp(27i7/3) (no sum on a). (9) 

The strings T a in (8) correspond to the numbered strings in Fig. 6(b). As 
in Section IIIA the triple string in this figure is a phantom string, with no 
Dirac string or thick vortex, so it contributes nothing to the action. 
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We will compute the action of the nexus, that is, the total action as- 
sociated with (8) minus the action associated with the vortices. Only the 
special case where the non-trivial part of string a runs along the positive 
coordinate axis a all the way to infinity will be given explicitly In eval- 
uating the action from (1) observe that the Tr\DiU\ 2 term is the same as 
Tr\Ai~ UdiU^ 1 ^ = Tr\Ai\ 2 , and that this form involves only the short-range 
Am part of the A part comes from UdiU^ 1 . This mass part of the action 
contains the short- distance logarithmic divergence mentioned earlier, but it 
occurs only in the thick vortices and not in the nexus, which has finite action. 

Taking the traces in the action (1) yields the action 



where 



and 



l(ab) 



1 

3? 



I = J (12) + 7 (13) + J (23 ) (10) 
d 3 x{[B^ - B^f + M 2 \A (a) - A {b) } 2 } (11) 



|W - M 2 f <— (#{*) qlq(a) ' 



1 ~ J (27r)2(gW-z e )(g2 + M 2) lC) g2 h V) 

A (a) _ f d3( l e^e ij{a) qi 
1 J (2vr) 2 (g( a )-k)(g 2 + M 2 )' 1 ' 

In these equations, q( a > and similar forms means the appropriate component 
along the axis a of the chosen string; e.g., q^ = q ■ e^, where e a is the unit 
vector along the string direction (assumed to be straight). In our case, with 
the string along the axes, one can drop the parentheses and interpret the a 
as ordinary spatial indices. The ie serves as a cutoff for the integration out 
to infinite distance along an axis, that is, one can think of the integration 
along an axis being stopped at a distance L = e _1 . The total action in (10) 
has a part diverging like e' 1 ; we will define the nexus action as this total 
action minus the divergent part. 

Straightforward calculation gives the total energy of the nexus as 

2vrM . , 

/ = -p-. (14) 



This is a remarkably small value, about one-tenth of the energy [17j of the 



QCD sphaleron. But unlike the sphaleron this quasi- Abelian nexus carries 
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no topological charge, as is easily checked. We will see, when we consider 
fully non-Abelian configurations next, that forming topological charge costs 
a considerable amount of energy. The energy of a nexus depends on the 
orientation of the strings, and the value in (14) is neither the smallest nor 
the largest energy as the orientations vary, but it is characteristic. Clearly, it 
is possible for the whole collection of vortices plux nexus to annihilate itself, 
by choosing the strings to coincide everywhere, but it is also possible to go to 
zero nexus energy in another way. Imagine taking strings 1 and 2 of Fig. 6b 
to lie along the negative z-axis, while string 3 lies along the positive z-axis 
as we started with. The configuration is geometrically like that of Figs. 1 or 
3 now. However, since Tn + T22 = — T33, one of the strings is reversed and 
this configuration turns out to be an ordinary center vortex with no nexus 
at all. 



IV. NON-ABELIAN NEXUSES 



The only interesting nexus in SU(2) gauge theory is not a quasi-Abelian 
one, which is just like the SU(3) nexus described above with strings 1 and 2 
coinciding. It is a truly non-Abelian object, roughly describable as a defor- 
mation of the SU(2) sphaleron mentioned above fl7f . There are also SU(3) 
non-Abelian counterparts. In both cases there is no analytic solution, and 
it would be tedious to find good numerical solutions, since one must solve 
several non-linear partial differential equations simultaneously (compare to 
the SU(2) sphaleron, which is described by a single non-linear ordinary dif- 
ferential equation). So we will content ourselves with some simple variational 
estimates for the nexus energy in SU(2); for SU(3) (and by a simple exten- 
sion SU(N)) we stop at setting up the (non-trivial) gauge kinematics for the 
nexus. The primary reason, aside from the complexity of going further, is 
that it is clear that the non-Abelian nexuses have considerably higher energy 
than the quasi-Abelian nexuses and are not likely to play an important role. 



A. SU(2) non-Abelian nexus 



Although the nexus has cylindrical symmetry only, it is useful to think of 
this nexus in relation to the usual Witten spherically-symmetric static ansatz 

iAj = -[($ - 1)X* + c.c] + f,J • rB i; (15) 
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Xj = e jak J a f k + i(Jj - TjJ ■ f). (16) 

Here J k = o"fc/2 are the £77(2) generators and $, Bi are functions of r. As 
usual $ = $!+2$ 2 is a complex scalar field coupled to a U (1) gauge potential 
B 1 whose r-component is B\. 

We can motivate further developments by referring to some well-known 
results on the Wu-Yang monopole, which is recovered at $, B\ — in (15). 
This monopole, because of a singularity at the origin, is not physically in- 
teresting, but it is mathematically simple. As is well-known this monopole 
can be transformed, by a singular gauge transformation, to an Abelian form 
with a Dirac string: 

A'j = UAjU- 1 + UdjU' 1 = l -J^j(l - cos6) (17) 

U = expiOJ-(j>; UJ-fU- 1 = J 3 . (18) 

This form is essentially the original Dirac monopole, with a string along the 
negative z-axis. The string flux is twice the minimum flux (corresponding to 
J = 2 in the notation of equation (3)), so it contributes nothing to Wilson 
loops. 

By a second singular gauge transformation with V = expz0J 3 we can 
shift half this flux to the positive z-axis (which replaces 1 — cos 9 by — cos 9 
in (17). At this point the magnetic fields associated with the vortex cores do 
point as desired (see Fig. 1). We can now, if we wish, make another gauge 
transformation with U~ 1 to restore spherical symmetry as much as possible. 

The combination of these three gauge transformations, namely 

U-i e i4>J 3U = e i4>J-f - w ^ ( 19 ) 

when applied to the Wu-Yang configuration, leads to a deformed kinematical 
ansatz which we give below in equation (21). Within this new ansatz we can 
describe the conventional QCD sphaleron [jlT , |I8| 1 and deformations which 
describe thick vortices emerging from it. 

The first step in implementing the new ansatz is to determine the appro- 
priate pure-gauge behavior of the nexus- vortex combination at large distances 
(by which we mean that both p and z are large). We compose this from the 
gauge W in (19) above, appropriate to represent the strings pointing into the 
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nexus of Fig. 1, and a pure-gauge form appropriate to a spherical configu- 
ration, such as a sphaleron. So we take the behavior at long distance of the 
gauge potential to be: 

Aj -> ZdjZ' 1 , Z = e mr)+<t>)J-^ (20) 

We now wish to apply these considerations to deform the usual QCD 
sphaleron JT7], This object is spherically-symmetric, has short-ranged 

field strengths, and a CS number of 1/2. It corresponds to choosing j3(r) in 
equation (20) to be 7T. The deformation we use is to be thought of as an ansatz 
for determining approximately the nexus mass by a variational principle; the 
ansatz is not necessarily one which fully expresses the kinematics of the 
nexus.0 So we write 

iAj = -[($ - 1)X* + c.c] + }jJ ■ fB x . (21) 

Note that in the last term on the right, the unit vector r has been replaced 
by the unit vector 0. Furthermore, the functions now depend not just on 
r, but on the cylindrical coordinates p, z, <ft, and the function B\ is taken to 
be 

B x = -- MKi(Mp), (22) 

that is, the vortex factor which multiplies in the center vortex taken by itself 
(equation (6)). Of course, this choice for B\ has cancellation of Dirac strings 
in its two parts, just as for the usual vortex. There is a difference with the 
vortex in equation (6), however; in (21,22) the SU(2) matrices enter as J -r, 
instead of J3 as in (6). This has the effect, as one readily checks, of giving 
Dirac-string singular field strengths proportional to J3 but with opposite 
sign on either side of the nexus, as schematized in Fig. 1. In the vortex, 
of course, the field strengths have the same sign of J3 wherever one goes 
on the vortex. Although this changes nothing about the asymptotic Wilson 
loop, which is indifferent to the sign of J3, it is a gauge-invariant distinction. 
There is no regular gauge transformation which changes J3 to J-f; the gauge 
transformation which does the job, in equation (18), is certainly singular. 

5 It is easy to write a cylindrically-symmetric ansatz which is kinematically self- 
consistent. This leads to coupled non- linear partial differential equations for nine scalar 
components, which we have not studied in detail. 
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The next step is to choose simple variational forms for the functions $j, 
which satisfy the appropriate boundary conditions at large distances from 
the nexus and vortex as expressed in equation (20), with /3 = tc, and lead to 
vanishing gauge potential at the origin. We will do this in analogy to such 
forms for the QCD sphaleron ||18| . The necessary dependence leads to: 



$1 = F(p,z) + G(p,z)cos0; (23) 
$ 2 = G(p,z) sin (j). 

The boundary conditions mentioned above lead to: 

p,z -> oo : F -> 0, G -> -1 (24) 
p,2 -> : F -> +1, G -> 0. 

Simple functions obeying these boundary conditions are: 

A^ + t 1 \ l + pr 

Here A is a variational parameter.^ The appearance of p in G is required to 
prevent singularities from derivatives of 0. 

It only remains to calculate the energy or action from equation (1), and 
to drop the vortex terms. The total energy has three terms: The first term 
depends only on $ 12 (equation (23)); the second term depends only on the 
vortex wave function B\ (equation (22)), and the third term is a cross term 
depending on all three functions. We drop the vortex term, and the resultant 
energy is the nexus energy E n . The G 2 part in equation (1) of the first term of 
the energy scales like 1/A, but the third (cross) term depends in a complicated 
way on both M and A if we use literally the expression (22) for Bi, where the 
mass scale is M. However, it turns out to make little difference if we replace 
M in (22) by the variational parameter 1/A, in which case the G 2 part of 
the cross term also scales like 1/A. For simplicity, we make this replacement. 
The nexus energy E n is defined as the total energy less the vortex energy, 
and after inserting equations (21,22,23,25) in the energy functional and doing 
the integrals we find: 

^ = 1.911^ + 2.197^^; (26) 

6 For the QCD sphaleron, the choice Bi—Q, $ = (A 2 — r 2 )(A 2 + r 2 ) gives a sphaleron 
mass only 1/2% larger than the true mass. 
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the M 2 X term comes from the mass term of the energy. Equation (26) has 
a minimum at A = 0.933/M, of value 3.22(4irM/g 2 ). This energy for the 
nexus is rather smaller than the conventional QCD sphaleron [|T7], which is 
about 5 .44(47r M / g 2 ) , and larger than that of the quasi- Abelian nexus, which 
is (see equation (14)) 2-nM/g 2 . 

B. SU(3) non-Abelian nexus 

We wish to find for SU(N) with N >2 the kind of nearly spherically- 
symmetric nexuses discussed above for SU(2). In particular, we desire to 
construct a nexus with N center vortices, each of unit flux, meeting at the 
nexus. It is clear from the properties of the Qi matrices (equation (4)) that 
this requires an embedding of SU(2) other than the standard ones, such as 
Ai/2, A2/2, A3/2 in S77(3); these standard embeddings cannot span the entire 
Cartan subalgebra, as the set of Qi does. Aside from the standard embedding 
of SU{2) as a subalgebra of SU(N), there is always another embedding, 
called the principal embedding f22|[] which does span the entire Lie algebra 
of SU(N). The kinematics of this sort of embedding have already been used 
in Ref. |16| for SU(3). Here we briefly review the principal embedding and 
the formation of a spherically-symmetric Cartan subalgebra, but stop short 
of a full-scale calculation of the properties of the resultant nexus. The reason 
for not carrying out the full calculation is that the rotation generators we 
must use for SU(3) are in the J = 1 representation of SU(2), and their 
trace is four times as large as the J = 1/2 generators used in the standard 
embedding of a unitary group in the rotation group SU(2). As a result, the 
nexus energy is about four times as large as the conventional QCD sphaleron 
energy, and is so large as to appear to be unimportant. Similarly, the CS 
number is four times as large, and has the value 2. 

In the principal embedding, the Lie algebra of the fundamental N x N 
representation of SU(N) can be formed from representations of SU(2) with 
spin J = (N — l)/2 as follows. Let J a be the generators of SU{2) with this 
spin, and form further elements of the Lie algebra of SU(N) via: 

M ab J a J b , M abc J a J b J c , ■ ■ ■ (27) 

where the M-tensors form a complete set of (numerical) symmetric and trace- 
less tensors and the series terminates at an M with N — 1 indices. Clearly 

7 I thank S. Ferrara for furnishing this reference. 
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the set (27) plus the J a form a Lie algebra with N 2 — 1 elements. One forms a 
spherically-symmetric version of the Cartan subalgebra by choosing a subset 
of the tensors M abc „, as the appropriate symmetric and traceless combinations 
of the unit vector f a , e.g., 

M ab = f a r b - U ab . (28) 

A spherically-symmetric Cartan subalgebra then is J • r, M a bf a fb, . . .. 

Now we specialize to SU(3). For this group the principal embedding 
is the familiar "nuclear physics" decomposition into 0(3) plus a quadrupole 
tensor. In this group one has { Ji, J2, J3} = {A7, — A5, A2} in terms of the Gell- 
Mann matrices and we normalize the symmetric and traceless quadrupole 
generators as: 

T ab EE \{J a , J b } - \l5 ah - ^Taa = 0. (29) 

The diagonal elements span the Cartan subalgebra. Explicitly: 

T n = diag(-2/3, 1/3, 1/3); T 22 = diag(l/3, -2/3, 1/3); 

T 33 = ^(l/3,l/3,-2/3). (30) 

The Lie algebra takes the form: 

[Ja, The] = itabdTdc + itacdTbd', (31) 
[Tab, Ted] = - J e ( ^ace^bd ^bce^ad ^ade $bc + tbde5ac)- (32) 

The T ab are normalized according to: 

1 2 

TrT ab T cd = -(S ac 5 bd + 5 hc 5 a d - -Sabred)- (33) 

The spherical basis for the Cartan subalgebra is J ■ r and 

Q = TahTab- (34) 

The most general spherically-symmetric form of the gauge potential Ai for 
SU (3) can be read off from the kinematic structure of the pure-gauge poten- 
tial 

Ai(r) = VdiV~\ V = expi(a(r)Q + f3(r)J ■ f). (35) 



16 



This yields: 

iAj = BxfjQ + dfjj ■ r + - 1)X* + c.c] (36) 

2 J 

where $ is a complex doublet scalar and Xj is the kinematic complex doublet 

Jj -fjJ • f + ie jM J k ri 
Kj + ie m K k ri 

and 

K i = rd i Q = 2{T ia f a -f i Q). (38) 

Here the functions Bi, C\, $ are all functions of r. This is, of course, the same 
basis as used by Corrigan et al. but in a different notation. It can be 
extended to time-dependent functions by adding A4 as a linear combination 

iA 4 = B 2 Q + C 2 J ■ f. (39) 

Under a gauge transformation given by V of (35), one finds that 

while _B 7 ,C 7 (7 = 1,2) are changed by gradients of as Abelian gauge 
fields would be. So the interpretation, similar to that of the corresponding 
Witten spherical ansatz for SU(2), is that of a complex d=2 scalar field with 
two i/(l) gauge symmetries as indicated by (40), with B, C the corresponding 
gauge potentials. 

Next we deform the spherically-symmetric ansatz as we did for SU(2), 
by constructing an SU(3) equivalent of equation (19). This requires first a 
specification of some angular functions whose gradient gives the long-range 
pure-gauge part of the vortices (the A term in equation (2)). The necessary 
generalization uses what is known, in the older electromagnetic literature, as 
the magnetic potential <3>. Corresponding to any closed curve T, the magnetic 
potential is defined as: 

If 1 

$r = ^ijk \ do-ijdk— — - (41) 
2 Js \x — z\ 

where the surface S has the curve V as its boundary. The gradient of $ gives 
the Ao term in (2), with a 2-k jump giving the Dirac-string singular magnetic 
field whenever the surface S is pierced by a closed loop linked to V . 
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We choose three curves, much as shown in Fig. 6(b), which has, as it 
must, two SU(3) nexuses joined by three thick vortices. There is also a line 
along which all the three strings 1,2, and 3 coincide, allowing for closure of 
these three curves. To construct a single SU(3) nexus in isolation, imagine 
that one of the two nexuses in the figure is taken to infinity. 

At large distances from the nexus and its attached vortices, the gauge 
potential must approach a pure gauge (c/. equation (20) for SU(2)). We 
choose this gauge to be: 

Ai -> RdiR-\ (42) 

R = exp*[Q($ 3 - - \§2) + \j ■ r($2 - $i)]. (43) 

One easily sees that along the line where the three strings coincide there is 
no 2tt jump in the exponent of R, and no Dirac string. On the other hand, 
the jump associated with any string taken by itself leads to a unit-flux Dirac 
string. For example, the 2n jump around string 1 leads to a flux factor 
expressed in 

2m 

R — > exp{-in[Q + J ■ f]} = exp — . (44) 

o 

If the strings are chosen to lie along the positive coordinate axes, as in the 
quasi- Abelian SU (3) vortex discussed above, the matrix R is unitarily equiv- 
alent to the same matrix with the coefficients of the $ a replaced by T aa ; this 
is the analog of the construction of the gauge matrix W for SU (2) in equation 
(19). 

V. SUMMARY AND CONCLUSIONS 

In SU (N) we have partially (fully, for the quasi- Abelian case) constructed 
some new solitons for the center-vortex picture, solitons which have up to N 
thick vortices extending from them, with total magnetic flux of zero (mod 
N). These solitons can be thought of as composed of N monopoles of unit 
flux. Such monopoles are unacceptable in isolation, as they would come with 
singular naked Dirac strings, but in the nexus these Dirac strings cancel. 
For N > 3 these solitons, if entropically favored, would lead to a vortex 
vacuum analogous to a branched-polymer gel. There is, however, no effect 
from this branching on the fundamental properties of the area law for Wilson 
loops. Some of the nexuses carry CS number (or topological charge in four 
dimensions). In fact, vortices themselves carry CS number or topological 
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charge, quantized in units of 1/N. Such quantization is quite natural from 



the point of view of solving the U(l) problem[23], but is not natural from 
the viewpoint of periodicity of the 9 angle in units of 2tt. However, these 
requirements can be reconciled, and we will discuss this in the center-vortex 
picture in a later publication. 

Even though the nexuses may not play a large dynamical role in QCD 
(at least at 9 = 0), it is worthwhile looking for them on the lattice, supplied 
with appropriate boundary conditions, as a test of the underlying properties 
of center vortices. It might be worthwhile to list here some of the other prop- 
erties of the center vortex picture which are accessible to lattice calculations 
and other tests. First, the baryonic Wilson- loop area law is predicted f24[ to 



be a A-law, not a Y-law. Second, the fundamental premise of the center- 
vortex picture is that confinement comes from a non-trivial group center. 
There are certain exceptional groups (the simplest is G2) which have only a 
trivial center, and it would be interesting to simulate them on the lattice to 
see if they confine. And third, there are tests yet to be devised which probe 
the interaction of vortices and the 9 angle. 
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FIGURE CAPTIONS 



Fig. 1. An isolated SU(2) nexus (circle) and unit-flux vortices (lines), 
with field strengths reversing at the nexus. 

Fig. 2. An isolated SU(3) nexus-vortex combination. 

Fig. 3. SU(2) nexuses in a closed vortex. 

Fig. 4. Possible configurations for SU(3) nexuses and vortices. 

Fig. 5. An isolated SU(N) vortex-nexus combination, with up to N 
vortices meeting at the nexus; the sum of fluxes must vanish (mod N). 

Fig. 6. Depiction of SU (3) vortex strings as closed Dirac strings. The 
double string in (a) is physical, but not in (b), since it carries no flux or 
energy. In (a) the double string is unstable to separation for the Abelian 
nexus. 
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